Domain Growth in Ising Systems with Quenched Disorder

We present results from extensive Monte Carlo (MC) simulations of domain growth in ferromagnets and binary mixtures with quenched disorder. These are modeled by the "random-bond Ising model" and the "dilute Ising model" with either nonconserved (Glauber) spin-flip kinetics or conserved (Kawasaki) spin-exchange kinetics. In all cases, our MC results are consistent with power-law growth with an exponent $\theta (T,\epsilon)$ which depends on the quench temperature $T$ and the disorder amplitude $\epsilon$. Such exponents arise naturally when the coarsening domains are trapped by energy barriers which grow logarithmically with the domain size. Our MC results show excellent agreement with the predicted dependence of $\theta (T,\epsilon)$.Comment: 11 pages, 15 figure

Kinetics of Phase Separation in Thin Films: Simulations for the Diffusive Case

We study the diffusion-driven kinetics of phase separation of a symmetric binary mixture (AB), confined in a thin-film geometry between two parallel walls. We consider cases where (a) both walls preferentially attract the same component (A), and (b) one wall attracts A and the other wall attracts B (with the same strength). We focus on the interplay of phase separation and wetting at the walls, which is referred to as {\it surface-directed spinodal decomposition} (SDSD). The formation of SDSD waves at the two surfaces, with wave-vectors oriented perpendicular to them, often results in a metastable layered state (also referred to as ``stratified morphology''). This state is reminiscent of the situation where the thin film is still in the one-phase region but the surfaces are completely wet, and hence coated with thick wetting layers. This metastable state decays by spinodal fluctuations and crosses over to an asymptotic growth regime characterized by the lateral coarsening of pancake-like domains. These pancakes may or may not be coated by precursors of wetting layers. We use Langevin simulations to study this crossover and the growth kinetics in the asymptotic coarsening regime.Comment: 39 pages, 19 figures, submitted to Phys.Rev.

Soliton Staircases and Standing Strain Waves in Confined Colloidal Crystals

We show by computer simulation of a two-dimensional crystal confined by corrugated walls that confinement can be used to impose a controllable mesoscopic superstructure of predominantly mechanical elastic character. Due to an interplay of the particle density of the system and the width D of the confining channel, "soliton staircases" can be created along both parallel confining boundaries, that give rise to standing strain waves in the entire crystal. The periodicity of these waves is of the same order as D. This mechanism should be useful for structure formation in the self-assembly of various nanoscopic materials.Comment: 22 pages, 5 figure

Chain length dependence of the polymer-solvent critical point parameters

We report grand canonical Monte Carlo simulations of the critical point properties of homopolymers within the Bond Fluctuation model. By employing Configurational Bias Monte Carlo methods, chain lengths of up to N=60 monomers could be studied. For each chain length investigated, the critical point parameters were determined by matching the ordering operator distribution function to its universal fixed-point Ising form. Histogram reweighting methods were employed to increase the efficiency of this procedure. The results indicate that the scaling of the critical temperature with chain length is relatively well described by Flory theory, i.e. \Theta-T_c\sim N^{-0.5}. The critical volume fraction, on the other hand, was found to scale like \phi_c\sim N^{-0.37}, in clear disagreement with the Flory theory prediction \phi_c\sim N^{-0.5}, but in good agreement with experiment. Measurements of the chain length dependence of the end-to-end distance indicate that the chains are not collapsed at the critical point.Comment: 13 Pages Revtex, 9 epsf embedded figs. gzipped tar file. To appear in J. Chem. Phy

Polymer Brushes in Cylindrical Pores: Simulation versus Scaling Theory

The structure of flexible polymers endgrafted in cylindrical pores of diameter D is studied as a function of chain length N and grafting density \sigma, assuming good solvent conditions. A phenomenological scaling theory, describing the variation of the linear dimensions of the chains with \sigma, is developed and tested by Molecular Dynamics simulations of a bead-spring model.Comment: 35 pages, 38 figure

Langevin Dynamics simulations of a 2-dimensional colloidal crystal under confinement and shear

Langevin Dynamics simulations are used to study the effect of shear on a two-dimensional colloidal crystal confined by structured parallel walls. When walls are sheared very slowly, only two or three crystalline layers next to the walls move along with them, while the inner layers of the crystal are only slightly tilted. At higher shear velocities, this inner part of the crystal breaks into several pieces with different orientations. The velocity profile across the slit is reminiscent of shear-banding in flowing soft materials, where liquid and solid regions coexist; the difference, however, is that in the latter case the solid regions are glassy while here they are crystalline. At even higher shear velocities, the effect of the shearing becomes smaller again. Also the effective temperature near the walls (deduced from the velocity distributions of the particles) decreases again when the wall velocity gets very large. When the walls are placed closer together, thereby introducing a misfit, a structure containing a soliton staircase arises in simulations without shear. Introducing shear increases the disorder in these systems until no solitons are visible any more. Instead, similar structures like in the case without misfit result. At high shear rates, configurations where the incommensurability of the crystalline structure is compensated by the creation of holes become relevant

Spontaneous creation of discrete breathers in Josephson arrays

We report on the experimental generation of discrete breather states (intrinsic localized modes) in frustrated Josephson arrays. Our experiments indicate the formation of discrete breathers during the transition from the static to the dynamic (whirling) system state, induced by a uniform external current. Moreover, spatially extended resonant states, driven by a uniform current, are observed to evolve into localized breather states. Experiments were performed on single Josephson plaquettes as well as open-ended Josephson ladders with 10 and 20 cells. We interpret the breather formation as the result of the penetration of vortices into the system.Comment: 5 pages, 5 figure

Transitions of tethered polymer chains: A simulation study with the bond fluctuation lattice model

A polymer chain tethered to a surface may be compact or extended, adsorbed or desorbed, depending on interactions with the surface and the surrounding solvent. This leads to a rich phase diagram with a variety of transitions. To investigate these transitions we have performed Monte Carlo simulations of a bond-fluctuation model with Wang-Landau and umbrella sampling algorithms in a two-dimensional state space. The simulations' density of states results have been evaluated for interaction parameters spanning the range from good to poor solvent conditions and from repulsive to strongly attractive surfaces. In this work, we describe the simulation method and present results for the overall phase behavior and for some of the transitions. For adsorption in good solvent, we compare with Metropolis Monte Carlo data for the same model and find good agreement between the results. For the collapse transition, which occurs when the solvent quality changes from good to poor, we consider two situations corresponding to three-dimensional (hard surface) and two-dimensional (very attractive surface) chain conformations, respectively. For the hard surface, we compare tethered chains with free chains and find very similar behavior for both types of chains. For the very attractive surface, we find the two-dimensional chain collapse to be a two-step transition with the same sequence of transitions that is observed for three-dimensional chains: a coil-globule transition that changes the overall chain size is followed by a local rearrangement of chain segments.Comment: 17 pages, 12 figures, to appear in J. Chem. Phy

Phase transitions in nanosystems caused by interface motion: The Ising bi-pyramid with competing surface fields

The phase behavior of a large but finite Ising ferromagnet in the presence of competing surface magnetic fields +/- H_s is studied by Monte Carlo simulations and by phenomenological theory. Specifically, the geometry of a double pyramid of height 2L is considered, such that the surface field is positive on the four upper triangular surfaces of the bi-pyramid and negative on the lower ones. It is shown that the total spontaneous magnetization vanishes (for L -> infinity) at the temperature T_f(H), related to the "filling transition" of a semi-infinite pyramid, which can be well below the critical temperature of the bulk. The discontinuous vanishing of the magnetization is accompanied by a susceptibility that diverges with a Curie-Weiss power law, when the transition is approached from either side. A Landau theory with size-dependent critical amplitudes is proposed to explain these observations, and confirmed by finite size scaling analysis of the simulation results. The extension of these results to other nanosystems (gas-liquid systems, binary mixtures, etc.) is briefly discussed